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1. Introduction and Outline 

Given the tremendous progress in high-precision cosmology, in particular, the decisive 
character of distance indicators and structure formation probes on large scales, the time 
is ripe to test the rigidity of Einstein's theory of gravitation on cosmological scales. This 
observational challenge is preceded by theoretical efforts aiming at consistent modifications 
of gravity at the largest observable distances. Obviously, only consistent theories are worthy 
to be confronted with data. 

In a classical theory, at the exact level, consistency refers to the existence of a well posed 
initial value formulation and continuous solutions for the underlying degrees of freedom on 
the entire spacetime manifold. More precisely, at the technical level, the evolution of a 
scalar degree of freedom $ on a spacetime manifold A4, should be given by a quasilinear, 
diagonal, second order hyperbolic equation 

(x; V$) V M V i/ $(x) = J (x; R; V$) , (1.1) 

where q is a smooth Lorentz metric, which, in general is not identical to the spacetime 
metric g, since it is permitted to depend on the scalar degree of freedom and its first 
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derivative, and J is a smooth function that may have a nonlinear dependence on these 
variables. Moreover, the current density J may depend on the Ricci tensor R(g). 

At the perturbative level, consistency of a classical theory demands hyperbolic evo- 
lution only on a bounded spacetime region, the perturbative domain, beyond which the 
fluctuation dynamics requires a non-perturbative completion that is consistent in the afore- 
mentioned sense. Perturbations around a classical solution can be quantized in the usual 
way, given technically natural interactions. The standard requirements for a probabilistic 
interpretation offer yet another and distinct notion of consistency related to the quantum 
stability of the theory. 

Classical stability at the perturbative level and quantum stability stand on quite dif- 
ferent footings. In fact, a finite domain of validity for classical perturbations does not 
cause any principal obstacle provided the underlying theory is consistent. Of course, once 
fluctuations leave the classical stability region their background develops an instability to- 
wards a new ground state. In contrast, a quantum mechanical instability is not related 
to specific initial conditions but instead to the massive production of particles at no cost, 
which are represented by negative norm states. Therefore, the underlying theory is flawed 
at the fundamental level. Additionally, what here is called quantum instability already 
has incisive effects within the framework of a purely classical analysis, which we discuss in 
Sec. 3.1. 

There are different frameworks for constructing consistent modifications of Einstein's 
theory of gravitation, once additional degrees of freedom are allowed in the description 1 . 
As an instructive example 2 , consider an additional second rank tensor not necessarily 
a metric, inducing the following relevant deformation of the Einstein-Hilbert action 



where H = g — m has mass dimension one and sets the characteristic scale for the 
deformation, and M(g) denotes the de Witt bimetric. Note that the de Witt bimetric 
is the most relevant albeit not unique choice for M, and we have neither written down 
explicitly the kinetic and potential self-interaction terms nor the matter sector. 

Assuming that VP is locked into the Minkowski metric, for one reason or another 3 , the 
interpretation of the deformation parameter follows from perturbing the metric around 
the Minkowski geometry, g = r/ + h. Expanding the action (1.2) to second order in the 
fluctuations h, the Fierz-Pauli theory [3] is rediscovered, for which the de Witt bimetric 
with respect to the background spacetime is the unique unitary choice. This justifies to 
think of the deformation as a mass term with the deformation parameter being the graviton 
mass. Of course, this interpretation hinges on the background geometry. 

The deformation presented in (1.2) was primarily investigated on Minkowski and de 
Sitter background geometries for the following reasons: Given the interpretation of the 

1 Hence, strictly speaking, these modifications are not faithful deformations in the BRST terminology. 
This is known as the statement that multi-diffeomorphic theories have no Yang-Mills analogue. 
2 For a bi-diffeomorphic construction see [1, 2]. 

3 For the moment it is not important to specify a dynamical mechanism that would give rise to the 
locking process. 
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deformation parameter on a Minkowski background, (1.2) has been used to study conse- 
quences of a graviton mass for the principle of equivalence, in particular, how the impact of 
seemingly technically unnatural sources on the background geometry could be weakened. 
Higuchi [4] showed that an intriguing relation between the deformation parameter and the 
cosmological constant needs to be fulfilled, m 2 > 1/3A = H 2 (H stands for the Hubble 
constant), in order to render the free dynamics of h on a de Sitter geometry unitary. If this 
bound is violated, unitarity violating negative norm states are introduced in the respective 
Hilbert space. 

Both backgrounds are special in that no source specifications based on radiation or 
matter fields are required. This is of course different for generic Friedman cosmologies for 
which the Hubble parameter varies in time and, thus, the right-hand side of Higuchi's bound 
generalized to sourced Friedman geometries can be expected to become time dependent. 
In particular, it seems that for any deformation parameter at early enough times unitarity 
violation is inevitable. The observation that the Hubble parameter's flow backwards in 
time seems to induce quantum instabilities is a serious challenge for the viability of the 
considered deformation. In fact, it is not clear whether the theory (1.2) makes sense at all. 

In a recent paper [5] we have already addressed the question of generalizing the Higuchi 
bound to generic Friedman spacetimes. This investigation relied on the usual Stiickelberg 
completion of h in conjunction with the Goldstone boson equivalence theorem [2]. We found 
that the theory (1.2) is, naively, subjected to two distinct bounds on Friedman cosmologies 
characterized by time dependent Hubble parameters. One of them, 

m 2 >H 2 + H, (1.3) 

enforces the absence of negative norm states (unitarity bound), whereas the second, 

m 2 > H 2 + H/3 , (1.4) 

describes the region where hyperbolic evolution of the fluctuations is guaranteed (stability 
bound). Beyond this region, hyperbolicity breaks down. But this is no principal problem, 
since the breakdown is triggered by a strong coupling regime that simply invalidates the 
perturbative approach, demanding for a nonlinear completion. Now, for all reasonable 
Friedman sources, H < 0. As an important consequence, the (classical) stability bound 
imposes a stronger requirement on the deformation parameter than the unitarity bound. 
For concreteness, we assign a value to the deformation parameter such that the stability 
bound is satisfied for times t > t*. Evolving backwards in time, the (classical) stability 
bound will eventually be invalidated since the deformation parameter is constant for the 
most relevant deformation (1.2). This signals the onset of the nonlinear regime. Thus, 
the would be unitarity bound lies beyond the perturbative domain and its derivation using 
perturbation theory cannot be trusted. In this precise sense, the theory is self-protected 
against unitarity violations, and moreover, there is an open window of opportunity for a 
consistent nonlinear completion. 

Even though the Goldstone boson equivalence theorem represents a powerful diagnostic 
tool that allows to extract the leading short-distance behavior (and, futhermore, many 
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interesting phenomena related to the most relevant deformation of the Einstein-Hilbert 
term can be understood by employing it, as for example the structure of the Fierz-Pauli 
mass term, the vDVZ discontinuity [6, 7] or the Vainshtein radius [8], see also [2]), it applies 
only in normal neighborhoods characterized by sub- Hubble distances <l/m. 

The main purpose of the present paper is to extend our consistency analysis to the 
intermediate and low energy regime. The prime framework to achieve this is a full-fledged 
cosmological perturbation theory for all degrees of freedom. As usual, the metric fluctua- 
tions are decomposed into irreducible SO(3) tensors in accordance with the isometries of 
Friedman geometries. Compared to the m = case, the equation of motion for the second 
rank SO (3) tensor modes is deformed only by an additional hard mass term. This is due 
to the fact that the degrees of freedom carried by the second rank SO (3) tensor are gauge 
invariant in the undeformed theory. The equations of motion for the first (vector) and 
zeroth (scalar) rank SO (3) tensors change considerably in the deformed theory. This is a 
testimony of the fact that the deformed theory (1.2) apparently has no gauge redundancy. 
It should be noted, however, that the deformed theory has an equal amount of constraints 
compared to the gauge freedom possessed by the undeformed theory (and in fact could be 
understood as the gauge fixed version of the Stuckelberg extended theory). 

The importance of these efforts is easily illustrated by the following results: From the 
SO (3) vector sector arises a stability criterion that cannot be recognized by employing the 
Goldstone boson equivalence theorem. This additional criterion signals the presence of a 
tachyonic instability whenever 

k p 2 hys + 3H + 2m 2 > (1.5) 

is not satisfied. Here, k p h ys = k/a(t) denotes the physical wavenumber. On sub-Hubble 
scales, this criterion is always fulfilled and, thus, the dynamics extracted by employing 
the Goldstone boson equivalence is not affected by the tachyonic instability in the vector 
sector. In fact, the equivalence theorem does not cover this sector at all, as it is subdomi- 
nant compared to the scalar sector. In order to preserve stability on super-Hubble scales, 
however, we find the new bound 

m 2 >-3/2H. (1.6) 

For any choice of the deformation parameter, this bound will be violated in the sufficiently 
early Universe, and, as a consequence, the vector modes will develop a tachyonic insta- 
bility, thereby triggering the transition to a new ground state. This result supports the 
self- protection mechanism found and analyzed in [5]. The vector sector, thus, plays an 
important part in the stability analysis, although it does not participate in the Goldstone 
boson equivalence. 

The cosmological perturbation theory of (1.2) reveals more insight into the stability 
dynamics, even in the scalar sector. Most importantly, the unitarity bound (1.3) seems 
at work on all scales and not just on extreme sub-Hubble scales. Isolating the scalar 
sector, this poses a potential threat for the self-protection mechanism, since it is a priori 
not clear whether a strong coupling regime self-protects the theory also on super-Hubble 
scales. We have, however, shown analytically that the scalar sector is protected against 
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unitarity violations for k = in the same sense as it was for sub-Hubble domains. To be 
more precise, we again find a stability violating region that occurs before the system enters 
the would-be unitarity violating region when evolved backwards in time. Compared to the 
sub-Hubble case, this region is simply shifted to larger values of the time t, so it seems 
reasonable to assume that there exists such a stability violating region for all values of k. 
This conjecture is also confirmed by a numerical analysis. Moreover, as we have discussed, 
we know that the vector sector will become unstable whenever (1.6) is violated, and thus 
contributes importantly to the self-protection of the system. 

2. The evolution of small fluctuations in the deformed theory 

The deformed equations of motion for the metric field g following from (1.2) are given by 

G^(g) - m 2 M^f{g) H a p = -SirMf 2 T^(g, X ) , (2.1) 

where again H = g — ^f. ^ is assumed to be locked into some reference metric, by 
one mechanism or another. T denotes the energy-momentum source, which depends on 
matter and radiation fields %, the metric field, and, in principle, an effective cosmological 
constant, as well. Any solution of the undeformed Einstein equations will be respected by 
the deformation, provided is locked into the appropriate tensor. 

The Bianchi identity of the undeformed theory together with energy-momentum con- 
servation of the source implies the following four exact constraints on the combination 
H = g — in the deformed theory, 

VH^ - V V H = . (2.2) 

Consider now metric perturbations h = g — 7 around a Friedman background 7 com- 
patible with T. Assume VP to be locked into the respective Friedman metric and to be 
inert to the extend that it can be considered a fixed reference metric. Then H = h and 
the equations of motion for small /i-fluctuations following from (2.1) are 

SR^irf, h) - m 2 ( V + hi 7/ , A = -SvrMp- 2 M^( 7 ) 6T a0 , (2.3) 

to linear order in h. Here, 5R and 5T are the linearized Ricci and energy-momentum 
tensors, respectively. To this order, the four constraints are given by 

V " V,( 7 )7^V = > ( 2 - 4 ) 

which looks like a gauge constraint, but in fact is not. 

The spatial isotropy and homogeneity of Friedman backgrounds allow us to decompose 
the metric fluctuation h into irreducible tensors with respect to these isometries, 

h 00 = -E , (2.5) 
h i0 = a [OiF + d} , (2.6) 
hij = a 2 [A5ij + didjB + dyCy + D {j ] . (2.7) 
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Here, E, F, A, and B denote S0(3) scalars, Gi and Cj are the components of a transverse 
SO(3) vectors (d a G a = 0,d b Cb = 0), and the D^ denote the components of a transverse- 
traceless rank-2 SO (3) tensor {d a D ab = and S ab D ab = 0). 

The appropriate source for a Friedman spacetime is the energy-momentum tensor of a 
perfect fluid. Its perturbations can be decomposed in the same spirit 

5T o = Sp-ph 00 , (2.8) 
5T i = -{p + p)5ui+p h 0i , (2.9) 
6Tij = ph i + a 2 5ij5p , ( 2 -10) 

where the normalization condition g(u,u) = —1 and the background equation u M = 5q 
have been used. The three-velocity field 5u will be decomposed in a gradient and a curl, 
Su a = d a 5u + Sua- 

Using the irreducible SO(3) tensors from (2.5-2.7), the constraint (2.4) can be decom- 
posed accordingly, 

-3i -B + (A/a 2 )aF + SHE - 3HA - HB = , (2.11) 
(aF) - 3H(aF)\ - dj [E + 2A] = , (2.12) 
-(aGj) + ACj - 3H(aGj) = , (2.13) 



where B = AB. The constraint (2.11) is obtained from the v = part of (2.4), (2.12) from 
its v = i part proportional to a gradient of a scalar, and (2.13) from its v = i part given 
by a transverse vector. 

Now, we have all ingredients to linearize Eq. (2.3) and to equate the rank-2, 1,0 SO(3) 
tensor contributions separately. 

2.1 Rank-2 contribution 

The rank-2 SO (3) tensor contribution results from the transverse-traceless part of the 
spatial-spatial components of (2.3), and is given by 

-Dij - 3HDij + (A/a 2 ) - m 2 Aj = . (2.14) 

It is worth mentioning that (2.14) reduces to its counterpart in the undeformed theory 
in the m — > limit. This is a manifestation of the fact that the constraint (2.4) cannot 
support transverse-traceless modes and, as a result, general relativity can be continuously 
recovered in this sector. Provided the deformation parameter is small, m 2 < H 2 , the 
deformation term in (2.14) will not change the dynamics very much. In particular, the 
frozen mode on super-Hubble scales, —A/a 2 <C H 2 , is still present like in the undeformed 
theory. 

Concerning stability, the equation of motion (2.14) always yields stable solutions, since 
the coefficients of both, the Dij and terms coincide with the sign of the coefficient in 
front of Dij. As a consequence, displacements will always be pulled back to the equilibrium 
position. 

In the following, we will always use the same symbol for both the real space and Fourier 
space amplitudes of any dynamical variable like D^ . 
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2.2 Rank-1 contribution 

The deformed equations of motion (2.3) contribute two equations in the SO(3) vector sector 
of the theory, one from equating the spatial-temporal components, the other from equating 
the spatial-spatial components. As we will see, it suffices to consider the spatial-temporal 
equation together with the constraint (2.13) and momentum conservation to solve the 
vector sector. The vector part of the spatial-temporal equation is given by 

16^M p " 2 (p + p) ou v /a = (A/a 2 - 2m 2 ) G - (A/a 2 )aC . (2.15) 

For convenience, let us define Gj = aGj. From the constraints (2.13), it then follows that 

AC = G + (3HG) . (2.16) 
Inserting this equation into (2.15) yields 

16irM p - 2 (p + p) 5u v = -G - (3Hg\ - 2m 2 G + (A/a 2 )G . (2.17) 



A solution for the divergence-free part or the three- velocity field <5u v can be obtained from 
the momentum conservation statement in the corresponding sector, which is given by 

(p + p)5u v< j +3H(p+ p)5u v = . (2.18) 

This shows that the quantity (p + p)5u v oc 1/a 3 decays and can therefore be neglected at 
late times. As a consequence, the equation of motion for G (2.17) is source- free at late 
times. 

Investigating the stability of (2.17), we see that the Hubble- friction enters with the 
correct sign, whereas the terms with no time derivatives on G need to satisfy 

(A/a 2 ) + 3H + 2m 2 G>0 (2.19) 

to give a stable solution for G. Surely, in certain kinematical regions and for particular val- 
ues of the deformation parameter, the bound (2.19) will be violated, and, as a consequence, 
a tachyonic instability will be generated. Indeed, for sufficiently early times, there will be 
such an instability for all three- momenta, provided that H increases faster than — A/a 2 for 
decreasing t. This is the case, for instance, during radiation and matter domination, but 
not for the epoch when the cosmological constant dominates. In the latter case, the vector 
modes are always stable for arbitrary three-momenta. 

On extreme super-Hubble scales, — A <C (aH) 2 , the system develops instabilities when- 
ever the bound 

m 2 >-3/2H. (2.20) 

is violated. This bound is a new result that has not been obtained in the previous work 
[5] based on the Goldstone boson equivalence. The bound (2.20) is instrumental for the 
self-defense of the theory against unitarity violations: Consider an equation of state of the 



-7- 



form p(p) = wp ,w =const. For w < 10/3, the bound (2.20) is even stronger than (1.4) 
and, furthermore, supports the self-protection mechanism described in [5]. 

Once the equation of motion (2.17) for G is solved, the constraint (2.13) allows to solve 
for C up to a spatially homogeneous contribution which, anyhow, does not contribute to 
the spatial-spatial components of the metric perturbation, since C enters only with spatial 
derivatives. This clearly shows that the vector sector contains exactly one independent 
divergence-free three-vector field, and, thus, is inhabited by two independent degrees of 
freedom. 

2.3 Rank-0 contribution 

Like in the undeformed theory, the scalar sector is the most intricate. It contains as 
geometric ingredients the scalars A, B, E as well as F, and from the source Sp, Sp, and 
Su. Not all of these variables are, however, independent. Indeed, assuming a source with 
equation of state p = p(p) allows to reduce the dynamics to a set of two coupled second- 
order differential equations for A and B = AB: 

A = -3(1 - w)HA + w (A/a 2 ) A - [2m 2 - 6w (H 2 - m 2 /2)] A + 

+wHB + 2w (H 2 - m 2 /2) B + 
+HE - m 2 E(A, B) , 

(2.21) 

B = -7HB -4(H 2 + m 2 /2) B + 

-Y2HA - 3 (A/a 2 ) A - \2B 2 A + 

+ (12H 2 - A/a 2 ) E(A,B) , (2.22) 
where E is expressed in terms of A and B, 
H+(2- 3w) H 2 - m 2 ] E(A, B) = 

-(w- 1/3) HA-(w- 1/3) (A/a 2 ) A - \h + (1 + 6w) H 2 - (2 + 3w) m 2 j A + 

- (w - 1/3) HB - (1/3) H + (1 + 6w) H 2 - (2 + 3w) m 2 B. (2.23) 

The remaining geometrical SO(3) scalar F can be obtained using the deformation constraint 
(2.11). Then 5p can be derived from the temporal-temporal component of the linearized 
deformed equations of motion (2.3), and Su can be derived from the spatial-temporal 
components of (2.3) by extracting the spatial gradient contributions. Finally, 5p follows 
from the equation of state 5p = c 2 5p where c s denotes the isentropic sound speed in the 
source. The details of this calculation can be found in the appendix. 



3. Stability analysis in the scalar sector 

In [5] we have already discussed some qualitative differences between the two bounds (1.3) 
and (1.4): The former leads to negative norm states, which spoils the probabilistic inter- 
pretation of the theory, while the latter signals the breakdown of perturbation theory. In 
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Sec. 3.1 we reiterate on the issue by presenting further arguments for the physical difference 
of both bounds, based purely on the classical evolution. After Sec. 3.1 we continue with 
the stability analysis in the scalar sector. 

3.1 Classical effects of the different types of instabilities 

As already mentioned above, what here is called quantum instability (that is the appearance 
of negative norm states in the quantized theory) already has an incisive effect within the 
framework of a purely classical analysis: Let us have a look at a setup, which is actually 
capable of capturing all the relevant physics at the linear level for sub-Hubble scales [5], 
based on the classical equation of motion for a scalar <fi, cup + e<j) + f34> = 0. Here, the 
coefficients a, f3, and e are functions of time. For a, (3, e > the system is stable. The 
classical stability bound manifests itself in a change of the sign of /3 while a is still positive, 
which triggers an exponential instability, and the perturbative analysis breaks down. For a 
gradual zero-crossing the spring constant is already small before the hard bound is hit and 
the oscillations might enter the nonlinear regime already before the exponential instability 
is triggered. 

Nevertheless, we can still choose initial conditions that allow us to evolve the system 
for a small amount of time inside the region f3 < until the fluctuation grows large. We 
can, however, not use this approach to try to cross the point where a turns negative as well, 
as close to this point, a is already small, and the effective spring constant has an extremely 
negative value, which goes to — oo just at the zero-crossing. Hence, in its vicinity, the 
time for which we can evolve the system in the just described fashion goes to zero. As 
a consequence, there is no reason why a change of sign of a after a change of sign of 
(3 should have any physical relevance for the full system. This is a manifestation of the 
self-protection mechanism. 

Let us now consider the opposite case when a changes its sign before f3 does. In this 
case, the effective spring constant f3/a grows big before the zero-crossing of a, confining the 
oscillations of 4> to small values even more. The equation of motion, however, runs into a 
singularity because the term with two time derivatives (thus terminating the time evolution 
of the system) vanishes. Hence, this case would be much more severe, as the system cannot 
even be evolved across the point where a vanishes. A possible counterargument to this 
reasoning is that the system enters the strong coupling regime whenever a — > 0. We will 
argue, however, that the described singular behavior of the equation of motion persists 
in the same way in the non-linear theory: Consider the non-linear term ^(j><f) that will 
become important once a ~ <^y. In fact, this is the only relevant non-linear contribution, 
since any other term containing two time derivatives but more fields, such as </> 2 <?i>, will be 
subdominant due to the fact that 4> itself is small, as explained. Thus the combination 
(a + 7</>)</> will determine the time evolution of the system, with the equation of motion 



Again, as long as (a + j(f>) > 0, the effective spring constant of the system grows large 
and confines (ft to small values. At best, 7^ might have some positive value, so that a can 



(a + j4>)(/) + e<f> + P<f> = 0, 
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become negative, but now a eventually drops to large negative values and will certainly 
overshoot the contribution 70 which is still small due to the small fluctuations. Hence, 
even the sum a + 70 will pass through zero and result in a singularity of the system. 

Let us elaborate a little bit more on the question why a vanishing coefficient a + 70 
in front of the term entails an unacceptable singularity. We will name the time of zero 
crossing to, that is 

a(t ) + 7(to)0(to) = 0. (3.2) 

Assuming that is regular at to yields the constraint e(to)0(to) + /3(to)0(^o) = by virtue 
of the equation of motion (3.1). Moreover, (3.2) yields the additional constraint 0(to) = 
— a(to)/7(to)- These constraints completely spoil the Cauchy problem as they allow only 
one particular choice of initial conditions. This clearly illustrates the singular behavior of 
(3.1) under the assumption of regular 0. 

Thus, we try to abandon the assumption of regularity of 0, and instead assume that 
~ (a + 70) _1 around to- Taylor expansion of the vanishing coefficient gives the leading 
behavior ~ (t — to) _<5 - The case 5 = 2 results in ~ ln(|t — to|) which is singular at t = to 
and thus unacceptable. The same is true for 5 > 2, for which we obtain ~ (t — to) _5+2 . 
If instead we have (5 = 1,0 would behave as ~ (t — to) hi(|t — to|) — (t — to), which would 
be well-defined at t = to- The term e0 in (3.1), however, would still be singular for this 
behavior of 0, such that this behavior cannot give a solution to the equation (3.1). 

3.2 Unitarity Bound 

At the level of the action for the SO (3) scalar A, the sign of the prefactor in front of the A 2 
term is crucial for the absence of negative norm states. (See [5] for details.) At the level of 
the equation of motion, this sign is determined by the prefactor of the A term which can 
be derived from combining equation (2.21) with the corresponding prefactor in the E term 
from (2.23). Combining both prefactors gives 

•" 2 - fl2 -" -A. (3.3) 



(l-a ! |l + 3c,') 

Evidently, in the scalar sector unitarity seems to require that m 2 > H 2 + H, which is 
precisely the bound (1.3) found in [5] by employing Goldstone boson equivalence. As an 
important result, we re-derived this unitarity bound in a full-fledged cosmological perturba- 
tion analysis, with a very important qualification: we find that the unitarity bound applies 
at all energies, and not just in the high-energy regime considered in [5]. 

In the following we solve the coupled equations of motions (2.22, 2.21) for the scalars 
A, B numerically, and analyze the stability of these solutions. For clarity, we subdivide the 
kinematical domain in three subdomains: extreme sub-Hubble scales (k 2 /a 2 3> m 2 ,H 2 ), 
intermediate scales, and extreme super-Hubble scales (k 2 /a 2 <C m 2 ,H 2 ). 

3.3 Extreme sub-Hubble scales 

This regime has been investigated previously [5] employing the Goldstone boson equivalence 
as a diagnostic tool to extract the leading short-distance dynamics. 
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From the full, coupled set of linear differential equations (2.21-2.23) these dynamics 
can be recovered by means of the adiabatic ansatz A, B,E oc e^, which is best for large 
kph ys : Introducing the ansatz into the system of equations and solving the (biquadratic) 
secular equation C4A + C2A 2 + c$ = 0, which results to leading order in large k p h ys , yields 

\/2A = ±yj (— C2 ± \J c\ — &CQCi)/ci, where A 2 = /i 2 |k p h ys | 2 - (Upper and lower signs can 
be chosen independently, which leads to four combinations.) In order to have a stable 
system, none of the eigenvalues may have a positive real part. Therefore, the presence of 
the outer ± implies that all eigenvalues must be purely imaginary. That necessitates that 
c 2 . < 4coC4 and that Co, C2, and C4 must have the same sign. Unitarity requires further 
that c 4 = 2{H 2 + H - m 2 ) is negative, reproducing Eq. (1.3). Hence, the system is stable 
when all coefficients are negative and c 2 < 4coC4. Then, from co = (H 2 + H/3 — m 2 )w we 
reproduce Eq. (1.4) for w > 0. 

For w < this relation would be exactly the other way round, implying that the 
system would never be stable. This phenomenon is known already from unmodified general 
relativity [9], where a system filled by a perfect fluid with w < is always unstable as long as 
kphys is not very small. As it is already present in general relativity, this instability cannot 
have anything to do with the degree of freedom used in the Goldstone boson equivalence 
analysis, which is absent in general relativity. This explains why said instability goes 
unnoticed in this case. It is important to notice that in this respect a scalar field does not 
correspond to a perfect fluid [10], which explains why this bound is also not obtained in 
[11]. 

Coming back to w > 0, the bound derived from C2 is always weaker than the stability 
bound (1.4), which follows from cq, or the requirement c 2 , < Acqc^. For w > 1/3 the 
requirement that c\ < ^cqc^ would be stronger than the bound (1.4). A numerical analysis 
in the regime where c\ < 4coC4 shows, however, that there is no instability as in the case 
where (1.4) is violated. While the latter leads to a clear exponential explosion forwards and 
backwards in time, the latter manifests itself in a beat with an amplitude of the envelope 
that grows relatively mildly backwards in time. Here the requirement c\ < 4coC4 obtained 
in the framework of the adiabatic analysis does not seem to give a relevant bound. Also in 
the case where the condition C2 < is violated, numerically no instability can be detected. 

Figure 1 shows the numerical solution for the scalars A and B in a radiation-dominated 
universe (c 2 = 1/3). The parameters were chosen such that (units unspecified) fc p h ys = 
250 /y/i, m = l/\/l~2, and H = l/(2t). Hence, k 2 /a 2 » m 2 ,H 2 is guaranteed for times 
t G [0.8,2]. The initial conditions have been chosen at t = 2, such that the system is 
evolved backwards in time. 

Let us first investigate the behavior of B. For times t > 1, B is oscillating with 
a Hubble-damped amplitude, clearly showing a healthy hyperbolic evolution forward in 
time. Evolving backward in time, however, B develops an instability for t < 1. Indeed, 
the parameters have been chosen such that the stability bound (1.4) is violated for t < 1 
This confirms the results of [5]. The behavior of A is similar, except that it develops the 
instability at an earlier cosmological time scale (which is later from the point of view of the 
system evolving backwards in time) , and oscillates with a higher frequency as compared to 
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High energy regime High energy regime 




0.8 1.0 1.2 1.4 1.6 1.8 2.0 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

Time t Time t 

Figure 1: Scalars A and B during radiation domination deep inside the Hubble radius, 
the scalar B. 

The basic properties of the solution are independent of the source's equation of state 
in the interval < cj 2 < 1. The case of a de Sitter source (cj 2 = —1) is borderline, since the 
parameter range for which the classical instability is triggered coincides precisely with the 
range of parameters for which unitarity gets violated. Hence, the strong coupling regime 
goes hand in hand with negative norm states. (See [5] for details.) 

3.4 Intermediate scales 

Figures 2 and 3 show the solutions for the scalars A and B during radiation domination from 
intermediate to extreme super-Hubble scales, that is, for different values of the comoving 
wavenumber k or, equivalently, for the physical wavenumber k/a(t) at time t = 1. For 
convenience and clarity, the other parameters have been chosen precisely as in the previous 
section. Like in the previous case, the initial conditions have been chosen at t = 2 and the 
scalar modes have been evolved backwards in time. For concreteness, the initial conditions 
are given by A = 0.01, B = 0.01, and dA/dt = 0, dB/dt = at t = 2. Note that the 
qualitative behavior of this dynamical system is quite insensitive to the choice of initial 
conditions, in particular, with respect to the stability analysis. 

It can be seen that the scalar modes' behavior on intermediate scales (and also on 
extreme super-Hubble scales, see next section) is very different from the dynamics in a 
normal neighborhood (see previous section). Compared to the latter case, the instability 
triggered at t = 1 becomes less and less pronounced with decreasing wavenumber. In order 
to appreciate this fact, notice the different ranges of mode amplitudes covered on the y-axes 
in figures 2 and 3 as compared to figure 1. In fact, scalar fluctuations on super-Hubble 
scales show a power law behavior which is triggered by the cosmological singularity (i.e. 
by the singular coefficients H oc 1/t etc), and which is clearly distinct from an instability 
triggered by a non-hyperbolic evolution. 

3.5 Extreme super-Hubble scales 

In order to elucidate further this result, let us analyze the stability of the scalar zero modes, 
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Figure 2: Numerical solution of A(t) during radiation domination for different values of 
kphys{t = !)• The other parameters have been chosen to be the same as for Figure 1. 

which can be performed analytically. The zero modes of A and B satisfy (2.21-2.23), 
A = -3(1 - w)HA - [2m 2 - 6w (H 2 - m 2 /2)} A + 

+wHB + 2w (H 2 - m 2 /2) B + 
+HE — m 2 E(A, B) , 

(3.4) 

B = -7HB -4(H 2 + m 2 /2) B + 
-V2HA - \2H 2 A + 

+12H 2 E{A, B) , - 13 - (3.5) 









Figure 3: Numerical solution for B(t) during radiation domination for different values of 
&phys(^ = !)■ The other parameters have been chosen to be the same as in Figure 1. 



where E is expressed in terms of A and B as follows, 
H + (2- 3w) H 2 - m 2 ] E{A, B) = 

- (to - 1/3) HA - \H + (1 + 6w) H 2 - (2 + 3w) m 2 



A + 



(w - 1/3) HB - (1/3) H + (1 + 6w) i? 2 - (2 + 3to) m 2 



(3.6) 



- 14 - 



Instabilities in the extreme sub- and super-Hubblc cases 

io 5 i i 1 — I — I I I I 1 1 — I — I I I I 



LO 4 




10~ J 2 5 10 1 2 5 1 2 

w+1 



Figure 4: Instabilities in the extreme sub- and super-Hubble cases. In the orange region 
(top, detached), the system is classically unstable for k p h ys = 0. The dark-blue region 
(bottom, left) depicts the region, where unitarity would be violated. In the green region 
(adjacent to the former), the system is classically unstable for large k p h ys . 

As a consequence, in this limit, the system of two coupled differential equations for A and 
B reduces to a single equation of motion for the linear combination S = A + B/3, 

[C 2 (w; t)8 t 2 + Ci(w; t)d t + C (w; t)]S = Q, (3.7) 

where the coefficients C2,i,o depend on the equation of state parameter w of the source and 
on time via the Friedman background evolution. Explicit expressions for these coefficients 
can be found in the appendix. 

A sufficient condition for hyperbolic evolution on the entire Friedman manifold and 
thus, for classical stability, is given by C\jCi > and C0/C2 > for all times, for a given 
source equation of state parameter w. We can analyze how these stability conditions depend 
on the parameter w and time t. The result is shown in Fig. 4, where the orange region 
corresponds to the classical instability region for the zero mode S, and inside the dark-blue 
region unitarity would be violated. Figure 4 shows that for a source with equation of state 
parameter w ~ 0.11, the zero mode's dynamics is always stable, confirming our explicit 
numerical result for a radiation dominated Friedman universe discussed in the previous 
section. For smaller values of w, when evolved backwards in time, the zero mode will 
always first enter the region of classical instability (orange), which signals the breakdown 
of perturbation theory. Evidently, it cannot enter the unitarity violating region (dark- 
blue), without passing through the strong coupling regime (orange). For large momenta, 
the area of classical instability moves downwards and comes to rest exactly on top of the 
area where unitarity would be violated, which, thus, still cannot be reached without first 
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crossing the former (green). Hence, in this sense the strong coupling regime self-protects 
the scalar zero mode from unitarity violation, as well. As a consequence, it is not clear 
at all whether the thus diagnosed unitarity violating region is of physical relevance, as it 
lies well outside the perturbative regime. We can turn this argument around and conclude 
that no inconsistency is present within the perturbative regime. 



4. Conclusion 

In summary, using cosmological perturbation theory, we have proven the consistency of the 
most relevant Einstein-Hilbert deformation in the perturbative regime. The deformation 
itself achieves consistency via a self-protection mechanism that pushes potential unitarity 
violations beyond the weak coupling regime. This confirms previous studies concerning the 
deformation's nontrivial stability dynamics, based on a Stuckelberg completion of the de- 
formation in conjunction with the Goldstone boson equivalence [5]. Most importantly, this 
work extends the self-protection mechanism to encompass the entire kinematical domain, 
ranging from sub- to super-Hubble scales. 

It would be interesting to study the proposed non-linear theories [22, 23, 24] with a rigid 
FRW background to see whether they non-linearly exhibit the self-protection mechanism. 

As discussed in great detail in [12], the self-protection phenomenon is a prime example 
for the recently conceived classicalization mechanism [13, 14, 15, 16, 17] and extends it 
further to free field dynamics on curved backgrounds. 



5. Appendix 



5.1 Derivation of the evolution equations in the scalar sector 

To start with, consider the part of the momentum conservation equation JV^T^j 
is built up from a derivative dj of a scalar variable: 



5p + d ((p + p) Su) + 3H (p + p)5u+^(p + p)E 







that 



(5.1) 



We will specialize to an equation of state of the simple form 5p = ^5 p. By doing so, 
we restrict ourselves to the case of a one-component system. The more complicated case 
of multi-component systems can be investigated, but one needs further special information 
about the system (for example the separate energy-momentum conservation of each com- 
ponent if they do not interchange energy and momentum). Further, using the Friedmann 
equations, one easily shows that 8irG (p + p) = —2H. The fluctuation 5u can be expressed 
in terms of metric variables using the iO-equations of (2.3), where one again extracts the 
contributions built from a derivative of scalar variables, 



8irG(p + p)d j 8u = d j —HE + A — m 2 (aF) 
Using this in Eq. (5.1), together with equation (2.12), one derives 

8irG^5p -HE - (zH 2 + 2H^j E + m 2 E + A + 3HA + 2m 2 A 



0. 



(5.2) 



(5.3) 
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Since the spatial divergence of the bracket in (5.3) vanishes identically, we know that the 
expression in the bracket is equal to some function of time alone. As we know from the 
basic equation (2.3) that h^ u = 0, T^ v = (which corresponds to A = 0, B = 0, E = 0, 
5p = 0, etc.) must be a solution, this function of time must be identically zero. Hence, we 
obtain 



8irG-^5p -HE- [3H 2 + 2Hj E + m z E + A + 3HA + 2m 2 A = 0. (5.4) 

Next, we will consider the r/-equations of (2.3) from which we extract the part of the 
form didjS with S a scalar. This gives 

didj E + A — a 2 B - 3aaB - 2m 2 a 2 B + 2aF + 4dF = 0. (5.5) 

Using (2.12) we can reexpress 

dj (2aF + 4aF^j = dj faaF) + 2aF) = = dj (-4aF -2E- 4A) . (5.6) 

Inserting this in (5.5) and taking the trace of the result gives 



-B - 3HB - 2m 2 B - AH— F - ^E - 3^ A = 0. (5.7) 

a or a z 



Finally, using (2.11) we obtain, 

—B - 7HB - AH 2 B - 2m 2 B - YIHA - 3^ A - \2H 2 A - ^E + \2H 2 E = 0. (5.8) 

This equation is the first of the two basic evolution equations in the scalar sector, see (2.22). 
The 00-equation of (2.3) gives 

- 3^ ) Bp = —\HE - ^E - 3 ( H 2 + h) E + \m 2 E ■ 



+\A + 3HA+lm 2 A + 
+ \B + HB + \m 2 B 

~M aAF )- ( 5 - 9 ) 

Using (2.11) one can eliminate F from (5.9), 

4vrG ( 1 + 3^ J 5p= +\HE - ^E + 3H 2 E +\m 2 E 



dp J 



§1 - 6HA - 3HA - 6H 2 A + \m 2 A 



— \B — 2HB — HB — 2H 2 B + \m 2 B. (5.10) 
The jfc-equations proportional to Sjk give 

-4vrG (\ - Sp = \HE + (3H 2 + H^j E — \m 2 E 

-\A+ ^A - 3HA - \m 2 A + 

— \HB - \m 2 B 

+H^F. (5.11) 
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Let us again eliminate F using Eq. (2.3), 



AttG 1 



dp 



5p = \HE + HE- \m 2 E 

-\A +^A + 3H 2 A + -\m 2 A + 
+ \HB + H 2 B - \m 2 B. 



(5.12) 



Inserting this expression for 5p into (5.4) results in a second independent evolution equation 
in the scalar sector (2.21). 

Equating (5.10) and (5.12) allows us to eliminate 5p 



dp 

1 



HE + -HE 
2 



m 2 B 



1 



-m z E 



1 



A 

2a 2 ' 



A + ^A + 3H 2 A - -m 2 A + -HB + H 2 B + 



1 + 3| 



-HE-3H 2 E+-m 2 E 



-A + ^A - 3HA + -m 2 A 
2 2a 2 2 



+ -HB - HB + -m 2 B 
2 2 



(5.13) 



Here, in addition we have used Eq. (5.8) to eliminate B. 

Our ultimate aim is to express E in terms of A and £>. In the first place, one might 
think that Eq. (5.8) does the job for every mode k p h ys , but the problem with this equation 
is that the resulting expression for E would contain B, so that whenever E appears one 
would get three time derivates on B. This is something we should, if possible, try to avoid 
for the sake of tractability, and indeed, this is possible. One way (among others) is first 
to derive an additional equation in A, B, and E by just using the constraints (2.11) and 
(2.12): 



—F 

a 



^A(aF) 



-3HE + 3HA + HB + 3A + B 



1 



A(aF) - 2H—F 



-5H—F 
a 



:E 



a* a a a* 

15H 2 E - 15H 2 A - 5H 2 B - 15HA - 5HB 

-2%A. 
a z 



2^A 



A 



\E + 



Using in addition Eq. (5.8) to eliminate B this can be cast into the form 

3A + 6HA - 4-4 + 3HA + 3H 2 A - 3HE - 3HE + 
-HB + HB- 2m 2 B + H 2 B - 3H 2 E = 0. 



(5.14) 



(5.15) 



As it happens, the ratio of the coefficients in front of E and A coincides for equations (5.13) 
and (5.15). Therefore, by appropriately adding both equations, one eliminates E and A at 
once, leaving an equation, which can be solved explicitly for E in terms of A and B and 
their first derivatives. This equation is given by Eq. (2.23). 
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5.2 The evolution equations for k = 

In the case k p h ys = the equations of motion (2.21, 2.22, 2.23) reduce to 

22 + 22 2 (1 + 6w) - ^ (2 + 3w)} (A +±b) + 22 (A+ ±b) (-1 + 3w) 



12H 2 E- 12H 2 [A+-B 



-22 -H 2 (2- 3w) + 



12H[A + -B 



B - 3HB - m 2 B = 



(5.16) 
(5.17) 



22 2 (4 - 6w) + m 2 (l + {A + ^Bj + 22(7 - 3w) (A + ^Bj + 



m 



+ H 2 (-7 + 3w)-2H + — )E - HE + [A + -B 



(5.18) 



i.e. the equation of motion for S = A + \B (5.18) decouples, which we will abbreviate by 



C 2 {t)S + d{t)S + C {t)S = 
with Cb(i), C\(t) and C2(t) given by 
Co = 



(5.19) 



22 2 (4 - 6w) + m 2 1 + -w 



+ 



+ 



H 2 (-7 + 3w) - 222 + 



m 



m 



-22 



-22 - H 2 (2-3w) + 



m 



-22 - 22222(2 - 3w) 
Ci = (7 - 3w)H + 



22 + H 2 (1 + 6w) - ~ (2 + 3w) 



22 + 2ffff(l + 6tu) 

2 

m 



-H - H 2 (2-3w) + 



m 



22 + 22^(1 + 6tu) - — (2 + 3w) 



(5.20) 



+22 22 2 (-7 + 3w) - 222 + 



m 



-1 + 3w) [—22 -H 2 {2- 3w) + 



m 



22 {H + H l (1 + 6u>) - — (2 + 3w) 



-27 -H 2 {2- 3w) + 



m 



21 



-22 <j if (-1 + 3w) (-27 - 22 2 (2 - 3w) + ^-)" 1 - 22 (-1 + 3w) (-22 - 22222 (2 - 3to)) 



22 2 (-l + 3u;) 



C 2 = 1 - - 

-22 - 22 2 (2 - 3u;) + 

Acknowledgments 



-H-H 2 + \ 



-22 - 22 2 (2 - 3w) + 



(5.21) 
(5.22) 
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